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Networks having the geometry and the connectivity of trees are considered as the spatial support 
of spatiotemporal dynamical processes. A tree is characterized by two parameters: its ramification 
and its depth. The local dynamics at the nodes of a tree is described by a nonlinear map, given 
rise to a coupled map lattice system. The coupling is expressed by a matrix whose eigenvectors 
constitute a basis on which spatial patterns on trees can be expressed by linear combination. The 
spectrum of eigenvalues of the coupling matrix exhibit a nonuniform distribution which manifest 
itself in the bifurcation structure of the spatially synchronized modes. These models may describe 
reaction-diffusion processes and several other phenomena occurring on heterogeneous media with 
hierarchical structure. 



PACS Numbers: 05.45. +b, 02.50.-r 



I. INTRODUCTION 

In general, media that support nonequilibrium pattern formation processes are nonuniform on some length scales. 
Often the nonuniformity arises from the intrinsic heterogeneous character of the medium, typical of pattern formation 
in many biological systems, or from random discontinuities or from clustering in the medium. It is well known that 
heterogeneities can have significant effects on the forms of spatial patterns; for example they can produce reberberators 
in excitable media and defects can serve as nucleation sites for domain growth processes. Recently, there has been 
much interest in the study of dynamical processes on nonuniform networks, such as fractal lattices small world 
networks f|, hierarchically interacting systems j||5), random systems ||, etc. An especially interesting class of 
nonuniform geometries are trees whose hierarchical structure and lack of translation symmetry can give rise to a 
number of distinctive features in their dynamical and spatial properties. Examples of phenomena where hierarchical 
networks appear include DLA clusters, capillarity, chemical reactions in porous media fit], turbulence ecological 
systems Q, interstellar cloud complexes ]Tot |, etc. Hierarchical structures have also been studied in neural nets, 
because of their exponentially large storage capacity | pd[ . Although many hierarchical structures found in nature 
have random ramifications, here we study the case of simple, deterministic tree-like lattices. This allows to focus 
on the changes induced in spatiotemporal patterns as a result of the hierarchical structure of the interactions in the 
system. 

In this article we consider discrete reaction-diffusion processes occurring on general trees. The spatiotemporal 
dynamics corresponds to a coupled map lattice system defined on the geometrical support of a tree. In Sec. II, the 
coupled map lattice models for the study of general trees are introduced. A general notation for the treatment of 
deterministic trees is also defined. The diffusion coupling among neighboring sites of the lattice is described by a 
matrix, which exhibits an ordered structure. In Sec. Ill, the spectrum of eigenvalues and eigenvectors of the coupling 
matrix. The eigenvectors form a complete basis on which all spatial patterns can be expressed as a linear combination. 
Explicit calculations of the eigenvalues and eigenvectors are provided. Sec. IV presents a study of the bifurcation 
structure and stability of spatially homogeneous, periodic patterns on trees for a local dynamics described by the 
logistic map. Specific features emerge as a consequence of the ramified character of the lattice. Finally, a discussion 
of the results is given in Sec. V. 

II. COUPLED MAP LATTICE MODEL 

Trees constitute a class of hierarchical networks which can be generated by a process of successive branching starting 
from an initial element. In this paper, it is assumed that the branching rule, or ramification number R, is the same 
through the network. At the initial level of the tree, which we call level 0, there is one element or cell which splits 
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into R branches connecting R daughters cells, which comprise the level of construction 1. Each cell then splits into 
R daughters, producing R 2 sites at level 2. This construction continues until some level L, which we call the depth 
of the tree. There are R l cells at the level of construction, or layer, I, where I = 0, 1, 2, . . . , L. Thus, each cell in the 
lattice has one parent and R daughters, except for the level cell, which has no parent, and for boundary cells at the 
level I = L, which have no daughters. The number of cells lying on the boundary is R L . The total number of cells on 
the tree, or the system size, is N = (R L+1 — 1)/(R — 1). 

A cell belonging to a level I is connected only to neighbor cells lying in adjacent levels on the tree, i.e., to its parent 
and daughters. We do not consider direct interactions among cells belonging to the same level of construction. With 
these prescriptions, a tree is completely characterized by two parameters: the ramification R and the depth L. We 
shall use the notation (R, L) to indicate a tree possessing those parameters. 

Each cell in the tree can be specified by a sequence of symbols (ol\oli . . . a/), where I is the level to which the cell 
belongs. A symbol otk can be take any value in a collection of R different digits forming an enumeration system in 
base R, which we denote by {ei, e%, . . . , en}. The level cell can be assigned a different symbol, say (ccq) = (eo)- 

A cell identified with the sequence {a\a^ . . . ai), with I > 0, is always connected to its parent which has the 
label (a.\a.2 ■ ■ . a;_i), where the sequence of the / — 1 symbols is the same as in the first (I — 1) symbols of the 
daughter cell (ai«2 ■ ■ ■ on). If I < L, the cell (a±a2 ■ ■ ■ on) is also connected to its R daughters which are labeled 
by (aiaa ■ ■ ■ one\), {a\a% . . . 0^62), . . . , (a±a2 ■ ■ ■ onen)] where the sequence of the first I symbols is the same as in the 
mother cell (o>\ot2 ■ ■ ■ on). Figure 1 shows a tree with ramification number R — 3 and depth L = 3, indicating the 
labels on the cells. 

A tree might be considered as the spatial support of spatiotemporal dynamical processes with either discrete or 
continuous time We focus on reaction-diffusion and pattern formation phenomena on trees. The dynamical systems 
considered here are defined by associating a nonlinear function with each cell of a given tree and coupling these 
functions through nearest-neighbor diffusion interaction. In this way, a coupled map lattice describing a reaction- 
diffusion dynamics on a tree with ramification R and depth L can be expressed as 

1. Level cell: 



x t+ i(e ) = f(x t (e Q )) + 7 



y^a?t(e») - Rx t (e ) 



(1) 



2. Level < I < L cells: 



x t +i (ai ...ai) = f(x t (ai . . . aj))+ 
" R 

+7 12 x t(osi . . . a t ei) + x t {ot\ . . . - {R + l)x t {ai ... a;) 



(2) 



3. Level / = L cells: 



x t+ i(a x ...a L ) = f(x t (ai . . . a L )) + 7 [x t (ai . . . a L -i) - x t (a x . . . a L )} ; 



(3) 



where Xt{a\ . . . ai) gives the state of the cell labeled by (a%, . . . , on) on the tree at discrete time i; f(xt) is a nonlinear 
function specifying the local dynamics; and 7 is a parameter that measures the strength of the coupling among 
neighboring cells and plays the role of a homogeneous diffusion constant. This type of coupling is usually called 
backward diffusive coupling and corresponds to a discrete version of the Laplacian in reaction-diffusion equations. 

The above coupled map lattice equations can be generalized to include other coupling schemes, non- uniform cou- 
pling, varying ramifications within the network, or continuous-time local dynamics. Different spatiotemporal phe- 
nomena can be investigated on tree-like structures by providing appropriate local dynamics and couplings. 

Equations (0)-@ can be written in vector form as 

Xt+i =f(x t )+ 7 Mx t . (4) 

The state vector x t possesses N components x t (j),j = 0, . . . , N — 1, corresponding to the states x t (ai . . . on) of the 
cells on a tree (R,L) labeled with the sequence of symbols {a\ ... a;). The matrix M expresses the coupling among 
the components {xt(j)}. 

The components of x t may be ordered as follows. The level cell is assigned the index j = 0. All other cells labeled 
by (a% . . . on) can be associated to a unique integer index j = 1, . . . , JV — 1, by the rule 
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(ai 



.on) 



R l - 1 
R- 1 



a kR 



l-k 



(5) 



fc=i 



As an example, consider the cell labeled by the sequence (21) on the tree of Fig. (1). In this case R = 3 and the 
cell belongs to the level / = 2. According to the rule of Eq.(||), its index is j = 11. Its parent is labeled by (2), and 
has the index j — 3; while its daughters, labeled by (210), (211) and (212), are assigned indexes j = 34,35 and 36, 
respectively. 

Given this notation, the j— component of the vector-valued function f(x<) is [f(xt)]j = f(xt(j))- For a tree 
characterized by (R, L), the elements of its corresponding coupling matrix M, denoted by M(i, j), (i,j = 0,l,...,N— 
1), are 



M(i,j)=M(j,i) 



-R 

-(R 

-1 

1 





if i = and j = 
if i = j and < i < 
if i = j and i > B 1 



R-l 



(6) 



iii^j and (j = int(^) or i = mt(^)) 
elsewhere, 



where int(g) means the integral part of q. The matrix M is a N x N real and symmetric matrix. It should be 
emphasized that M plays the role of a spatially discrete diffusion operator on treelike networks, analogous to the 
Laplacian in a spatially continuous reaction-diffusion equation. 



III. SPECTRUM OF THE COUPLING MATRIX 

Similarly to reaction-diffusion processes on fractal lattices pL the spatial patterns that can arise on trees are 
determined by the eigenmodes of the coupling matrix M. Additionally, the stability of the synchronized states is 
related to the set of eigenvalues of M. 

In order to analyze the eigenvector problem, consider a tree with ramification R and depth L on which a spatiotem- 
poral dynamics has been defined in the vector form of Eq.(Q). The complete set of orthonormal eigenvectors of the 
corresponding matrix M can be described as the superposition of two distinct subsets of eigenmodes. One subset, 
which will be denoted by {u„}, comprises those eigenvectors associated to non-degenerate eigenvalues; and the other 
subset contains the eigenvectors corresponding to degenerate eigenvalues of M, and will be represented by {v^} (the 
indices refer to the degeneracy, as explained bellow). Thus, the complete set of eigenvectors of M is {u„} U {v^ ls }. 
Each eigenvector describes a basic spatial pattern that may take place on a tree characterized by (R, L). 



Non-degenerate eigenmodes 

The eigenvectors of M belonging to the non-degenerate subset {u„} satisfy 

Mu„ = fr„u„, n= 1,2,...,!/; (7) 

where b n is the eigenvalue associated to the eigenvector u„. There are v distinct eigenvectors with their corresponding 
eigenvalues in this subset. The j-component of a vector u„ is [u n ]j = u„(ai, Qf2, • ■ ■ , on), according to the associating 
index rule, Eq.(^). Any eigenvector u„ in this subset is characterized by the following property: all its components 
corresponding to cells of the tree lying on the same level or layer I are identical, i.e., 

u n (a>ia t ...ai) =u n ((3iP 2 ---f3i), (8) 

where (aa2 ■ ■ - cui) and (/3i/?2 ■ ■ ■ Pi) label any two cells belonging to the level I. A non-degenerate eigenvector of 
M thus possesses homogeneous layers. Because of property (ra), we also refer to the elements of {u„} as layered 
eigenvectors. 

An eigenvector u„ representing a tree of depth L has (L + 1) levels, including the level cell. Because of the 
homogeneous layer property, there can be v — L + 1 distinct eigenvectors in the subset {u„} satisfying this condition; 
one eigenvector for each homogeneous layer that can be formed. Note that the number of different layered eigenvectors 
depends only on the depth L of the tree and not on its ramification R. In particular, the spatially homogeneous 
eigenmode of M, which we denote by ui, belongs to the subset {u n } and its N components are 
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ui(aio 2 ...ai)=JV- 1 / a ; VI, Va fc (9) 

and 

ui = -j=col(l,l,...,l). (10) 
V AI 



Since the eigenvectors of M are mutually orthogonal, all other eigenmodes, in either subset {u„} or {v^ s }, must 
satisfy 

J2 u n (ai . ..ai) = 0; n^l, VZ, Vajt; 
E < s (ai...a ; ) = 0; VZ, Va fc ; 

ai, ...,aj 

that is, the sum of the components of any other eigenvector of M different that Ui is zero. 

Figure 2(a) shows the three layered eigenvectors and their associated eigenvalues of the coupling matrix correspond- 
ing to a tree with R — 3 and L = 2. 

The set of eigenvalues arising from Eq. (Q) may be ordered by decreasing value as &2, ...,£>£, bz+j.}- By 
Gershgorin's theorem [Q, the homogeneous eigenvector possesses the largest eigenvalue of M, which is b\ = 0. On 
the other hand, for large L the smallest eigenvalue is found to be 



-(2 + 3R)-y/R(A + R) 
hm b L+1 = — v - , (12) 

L — >oo Z 



and similarly, we find 



hm b2 =-(K + l) + VW±V^ 

L — >oo 2 

The eigenvalues 62, ... , 6l, &l+i} appear in pairs b n and & n ', as 6l+i and 62 above, according to the sign of the 
square root term. These pairs are related by 

b n + b n . = -2{R + 1), (14) 

where 

, / / £±2 , if i is odd 
[ L + 3 , if L is even. 

The eigenvalue 6(l+3)/2 arises whenever L is odd, and it is not associated with another b n . Its value is 

bL±3=-(R+l). (15) 
Thus, because of ( |l4| ) and (|l5|) the eigenvalues associated to the non-degenerate eigenvectors satisfy 

V 

J2K = -L(R + i)- (16) 

n=l 

Figure 3 shows the spectrum of eigenvalues {b n }, indicated by black dots, for a tree with ramification R — 3 at 
successive depths L. Eigenvalues associated to degenerate eigenvectors of M, to be discussed next, are also shown in 
Fig. 3. 

Degenerate eigenmodes 

The subset of degenerate eigenvectors {v^ ls } of the coupling matrix M corresponding to a tree characterized by 
(i?, L) satisfy 

MvL = a ms v B ms (17) 
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where a ms is the eigenvalue associated to a group of 8 degenerate eigenvectors {v,^ s , v^ s , . . . , v^ s } belonging to 
{ v msl- The index g goes from 1 to 8 and counts the different eigenvectors associated to the degenerate eigenvalue 
a rns . The integer indices m and s label different eigenvalues a ms . 

The j— component of a vector {v^ s } corresponds to a cell of the tree labeled by the rule i.e., [v^, s ]j = 
v ms( a ii a 2i ■ • • ) a i)- The eigenmodes {v^ s } are characterized by the following two properties, 

v° ns (aia 2 ...ai)=0; 1 = 0,1,2,..., L- m; (18) 

that is, all the components of v^ s lying in successive levels vanish up to the level I = L — M; and 

J2 vff ls (a 1 a 2 ...a i ) = 0; (19) 

at±ae2...a:i 

i.e., the sum of the components of an eigenvector v^ s , lying on the same level of the tree spatially described by v^ s , 
is zero. 

An eigenvector v^ s is spatially uniform in part, having all its components, or equivalent cells, equal to zero up to 
level L — m. The index m counts the number of remaining non- vanishing layers in the eigenvector vf ns , and its possible 
values are m = 1, 2, . . . , L. Each of the m non-vanishing layers may be homogeneous, but different among each other. 
The index s counts the number of possible different eigenvectors with m different homogeneous, non-vanishing last 
layers. Thus, s may take the values s = 1, 2, . . . , m. 

The index g lifts the degeneracy of vectors with the same indices m and s. The remaining, non- vanishing last layers 
may in fact be non-homogeneous, and may consists of subtrees with homogeneous sub-levels, which would reproduce 
the structure of the layered non-degenerate eigenvectors u„. The level I = L — m is the last vanishing layer in an 
eigenvector having m non-vanishing layers. On this layer, there are R L ~ m components or cells, and each of these cells 
gives origin to R layered subtrees, i.e., subtrees with homogeneous layers. These R subtrees themselves are related 
by the sum property Eq. (|l9|), which results in (R — 1) linearly independent subtrees. Therefore, the number of 
linearly independent eigenvectors v*^ s with the same values m and s is 8 = (R — l)R L ~ rn . The index g expresses the 
degeneracy of the eigenvectors associated to the eigenvalue a ms , and it may take the values g = 1, 2, . . . , (R — l)R L ~ m . 
In this way, the set of degenerate eigenvectors {v^ s } of the matrix M corresponding to a tree (R, L) is fully described. 

Figure 2(b) shows the subset of degenerate eigenvectors {v^ s } and the eigenvalues corresponding to a tree charac- 
terized by {R,L) = (3,2). 

The index m also expresses the form in which an eigenvalue a ms arises in the spectrum of eigenvalues of M. An 
eigenvalue a ms appears for the first time at a level I — m > 1 and stays in the spectrum of eigenvalues of M at 
subsequent levels up to / = L. Thus the index m may take the values m = 1, 2, . . . , L. On the other hand, m different 
eigenvalues arise at the level of construction I = m of the tree, which are counted by the index s = 1, 2, . . . , m. 

The total number of distinct eigenvalues of type a ms belonging to the spectrum of a matrix M associated to a tree 
(R, L) is 

L(L + 1) 

E m = 2 ■ (2°) 

rn—l 

Additionally, the m eigenvalues a ms that appear at a level / = m satisfy the following property 



Ea TOS = -(m-l)0R + l)-l, (21) 



and therefore the total sum of eigenvalues a ms for a matrix M associated to a tree (R, L) will be, 

L m L 



m—l s—1 rn — 1 



= _(fl + l)(L-l)L 

2 V ' 

Figure 3 shows the spectrum of eigenvalues {a ms } for a tree with ramification R = 3 at successive depths L. 
Eigenvalues {b n } corresponding to non-degenerate eigenvectors are also shown there. Thus the distribution of the full 
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spectrum of eigenvalues of the coupling matrix M can be seen as a function of the depth of the tree in Fig. 3. 
that the full spectrum {b n } U {a ms } is always contained between the eigenvalues b\ = and 

The total number of distinct eigenvalues of M, including both types a ms and b n , and denoted by f2, is 



0= (L + l) 



L(L + 1) (L + 1)(L + 2) 



Note 



(23) 



Note that the total number of eigenvalues in the spectrum of M is determined only by the depth of the tree and is 
independent of its ramification, although the specific values of the eigenvalues do depend on both R and L. Since there 
appear m eigenvalues of type a ms at each level / = m and there are (R—l)R L ~ m degenerate eigenvectors v^ s associated 
to each eigenvalue a ms , the total number of independent eigenmodes in the subset {v^,.} is J2m=i m (R — l)R L ~ m . In 
the non-degenerated subset {u n } there are (L + 1) independent eigenvectors, as we saw before. Therefore, the total 
number of independent eigenmodes of M is 



(L + l) + m (R- 



L-r 



(L + 1) + (R-1)R 1 



R 



L + l 



-R(L + 1) + L 



R L (R-1) 



R L+1 - 1 
R-l 



(24) 



as expected. 

Figure 4(a) shows the complete spectrum of eigenvalues of M and the degeneracy fraction of each eigenvalue, for 
a tree characterized by parameters (R,L) = (2, 11). The L + 1 = 12 eigenvalues b n are indicated by dots and they 
are non-degenerate, while the degeneracy S — 2 11_m of each of the L(L + l)/2 = 66 eigenvalues a ms is plotted as a 
vertical bar. It is evident that both the distribution of eigenvalues and their degeneracies are nonuniform. Another 
convenient representation of the scaling properties of the spectrum of eigenvalues of the coupling matrix can be 
obtained by plotting the accumulated sum of the degeneracies of all eigenvalues, that is the measure of the spectrum 
of M (denoted by p), on the eigenvalue axis for large L, as in Fig. 4(b). The resulting graph presents the features of 
a devils staircase, a fractal curve arising in a variety of nonlinear phenomena. 

The eigenmodes of the coupling matrix reflect the topology of the tree and they are analogous to the Fourier 
eigenmodes appearing in regular Euclidean lattices. In this sense, conditions Eq. (||) and Eqs. (18)-(p"9|) represent 
different wavelengths on a tree characterized by parameters (R, L). 



An example 



As an example of calculation of eigenvectors, consider any tree with ramification R and depth L = 1. 
of cells of the tree is N = (R 2 — 1)/(R — 1) = R+l. Thus, the tree consists of a mother cell at level I - 
to its R daughters at level I = 1. The corresponding (R+l) x (R. + 1) coupling matrix has the form 



M = 



The number 
connected 



-R 


1 


1 


1 . . 


• 1 \ 


1 


-1 





.. 


. 


1 





-1 


.. 


. 


1 








-1 .. 


. 



V i o o o 



(25) 



-1/ 



and the associated eigenvectors of M have (R + 1) components. There exist two non-degenerate eigenvectors, which 
arc the homogeneous Ui = col (1,1, . . . , 1), and U2, associated to the eigenvalues bi — and 62 , respectively. 

There is only the eigenvalue an associated to 5 = (R — 1) degenerate eigenvectors {vj 1; v^ 1; . . . jvfj -1 }. The total 
number of independent eigenmodes of M is (R+l), and the total number of distinct eigenvalues is f2 = 3, in agreement 
with Eq.ph. 
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All the eigenvectors of M have the level 1 = 0, with one component or cell; and the level I = 1, for which there 
are R components. The layered eigenvector u 2 will have the form 112 = col(x, y, y, . . . , y). Its eigenvalue equation 
MU2 = b 2 u 2 , plus the normalization condition |u2| = 1, yield the relations 



-Rx + Ry = 622;, 
x-y = b 2 y, 
x 2 + Ry 2 = 1, 



(26) 



whose solutions are b 2 = —(R+ 1), x 



, - R y= , 1 
1 



Thus, 



U2 



y/WR + T) 



col (-R,l, 1, •••,!)• 



(27) 



On the other hand, the degenerate eigenmode has the form = col (0, x±, x%, . . . ,xr), satisfying properties 
( |l8| ) and ([l9|), as well as the eigenvector equation Mv'j = auvfi and the normalization condition jvf-J = 1. These 
relations lead to 



1, 



T R x 2 



(28) 



which imply that an = —1. Making x\ = x and x 2 = X3 



xr = y, we get 



x+(R-l)y = 0, 
x 2 + (R - l)y 2 = 1, 



(29) 



with solutions, x = ,[ R _ x \ , y = 



Thus 



'11 



: COl 



0,-(i?-l), 
V (-R-i) times/ 



(30) 



For = col (0, 0, X2, ■ ■ ■ , xr) the procedure can be repeated by making x% = Z, £3 = X4, = . . . = xr = w, obtaining 



'11 



col 



0,0,2:, w, . . . , w 
\ (i?.-2)timesy 



(31) 



z+(R-2)w = 0, 
z 2 + (R-2)w 2 = 1, 



In general, we get 



-OR -2) 



y/(R-l)(R-2Y 



y/(R-l)(R-2) 



vf x = col 



0, ...,0,2, 

\g times (R-g) times/ 



(32) 



z + (R - g)w = 0, 
z 2 + (R-g)w 2 = 1, 
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z = — ( R 9 2 w= - 1 g = i,2,.. -,(-R-i); 

giving (i? — 1) eigenvectors in the degenerate subset {v^ s } for this example. With the addition of the two non- 
degenerate eigenvectors Ui and U2, there are N = R + 1 independent eigenvectors. Therefore, all the eigenvectors 
and eigenvalues associated to a matrix M corresponding to a tree characterized by parameters (R, 1) are accounted 
for. Note that since the eigenvalue an = —1 appears at level I = m = 1, it will stay in the spectrum of eigenvalues 
of M for all subsequent levels of construction, i.e., an = —1 arises for trees of any depth. Similarly, the eigenvalue 
bi = and its associated homogeneous eigenvector Ui always appear in a tree. 



IV. BIFURCATION STRUCTURE AND STABILITY OF SPATIALLY SYNCHRONIZED STATES 

Spatially synchronized states in extended systems are relevant since we are often interested in the mechanism by 
which a uniform system breaks its symmetry to form a spatial pattern as a parameter is changed. Consider spatially 
synchronized, period K states such as Xt(a\, . . . , a;) = Xk, V(ai, . . . , a;); where Xk ,{k = 1, 2, . . . , K) is a period K 
orbit of the the local map, f^ K \xk) — %k- The linear stability analysis of periodic, synchronized states in coupled 
map lattices is carried out by the diagonalization of M in Eq. (Q) , and it leads to the conditions |l3| 

K 

II [/'(**) +7M]=±1, (33) 

k=l 

where /i is an eigenvalue, in either set {b n } or {a ms }, of the coupling matrix M describing a tree (R,L). There are 
Cl=(L+l)(L+ 2)/2 different values of ft (Eq. @) to be used in Eqs. @. 

The nonuniform distribution of the eigenvalue spectrum is manifested in the stability of the synchronized states 
through this relation and give rise to important differences when compared, for instance, with the bifurcation structure 
on regular lattices. As an application, consider a local dynamics described by the logistic map, f(x) — Xx(l — x). In 



this case, the bifurcation conditions, Eq. (33), for the period K = 2 P , synchronized state on a tree characterized by 



parameters (i?, L) can be expressed as the set of curves 

v 

SIM = J] [A(l - 2x k ) + 7 /i] - ±1- (34) 
fe=i 

For each sign, Eqs. ( |34| ) yield (L + 1)(L + 2)/2 boundary curves in the plane (7, A) which determine the stability 
regions of the period 2 P , synchronized states on the tree. 

The scaling structure for the period-2 p , synchronized states in trees is similar to that of a any lattice described by 
a diffusive coupling matrix, since the form of Eq. (|34|) is the same in any case. As for any lattice (for example regular 
Euclidean lattices Q or fractal lattices 0]), the stability regions for the period-2 p , synchronized states in the (7, A) 
plane scale as A ~ 6~ p , and 7 ~ a~ p , where 6 — 4.669 . . . and a = —2.502 . . . are Feigenbaum's scaling constants for 
the period doubling transition to chaos. However, the specific structure of the eigenvalue spectrum of the coupling 
matrix determines the shapes and gaps of the regions of stability of synchronized, periodic states. 

The boundary curves Eqs. ( |34] ) for the synchronized, fixed point state (j> = 0) are given by the straight lines 

A = /i7 + l, A = //7 + 3; (35) 

which are first crossed for the most negative eigenvalue, fj, = Fig. 5(a) shows the boundary curves Eq. (^) in 

the plane (7, A) for the period-two (p = 1), synchronized state on a tree with (i?, L) = (3, 3), which are given by the 
two sets 

Sl(b n ) = -A 2 + 2A + 4 + 7 6„( 7 6„-2) = ±1; 71 = 1,2,3,4; (36) 

and 

Sl(a ms ) = -A 2 + 2A + 4 + 7a ms ( 7 a ms - 2) = ±1; (37) 
to = 1, 2, 3; s — 1, . . . , to. 

The boundary between the synchronized, fixed point state and the synchronized period- two state is at A — 3. The 
upper boundaries (corresponding to —1 in the r.h.s of Eqs. (p6])-(p7|)) have minima A m i n = 1 + \/h at values 7 m i n = 
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l/b n and 7 m i n = l/a ms (for any period 2 P , A m ; n depends on p). Fig. 5(b) shows a magnification of Fig. 5(a) 
around the minima of the upper boundaries. The distribution of the minima 7 m i n and the presence of nonuniformly 
distributed gaps (niches) in the boundary curves reflect the nonuniform structure of the eigenvalue spectrum. Since 
the nonuniformity in the distribution of eigenvalues persists at any depth L of a tree, this property allows for regions 
of stability of the synchronized states (niches) characteristic of trees and which are not present in other geometries, 
for example in regular lattices, where the distribution of eigenvalues of the coupling matrix is uniform and continuous 
in the limit of infinite size lattices. 

The set of eigenvectors {u„} U {v^ s } of the coupling matrix M constitute a complete basis (normal modes) on 
which a state x t of the system can be represented as a linear combination of these vectors. The evolution of x t then 
reflects the stabilities of the normal modes. Fig. 5(b) shows how the synchronized state may become unstable through 
crossing of the upper boundary; the first boundary segment crossed determines the character of the instability. For 
example, consider an initial state consisting of a small perturbation of the synchronized, period-2 state at parameter 
values just beyond the boundary segment corresponding to a.22 indicated by a cross in Fig. 5(b), where this initial 
state is unstable. The inhomogeneous period-4 final spatial pattern is represented in Fig. 6; it corresponds to a linear 
combination of the six cigenmodes {vf 2 ; g = 1, . . . , 6} and associated to the degenerate eigenvalue (Z22 of the matrix 
M corresponding to the tree (R, L) = (3, 3). All other modes are unstable in this region of parameter space. For any 
depth L of the tree, and any period 2 P , the boundary curve 63(0122) = — 1 separates a niche of the synchronized state 
from the stable region for the eigenmodes vf 2 corresponding to a 2 2- Thus, a transition between these two spatial 
patterns can always be observed in the appropriate regions of the (7, A) plane. 

V. CONCLUSIONS 

In a system of interacting agents, such as the models presented in this article, the coupling matrix contains the 
connectivity of the network and it determines the spatial patterns that can arise in the system. The underlying 
inhomogeneous structure of trees has pronounced effects on the spatial patterns that can be formed by reaction- 
diffusion processes on these lattices. The spatial patterns are determined by the eigenvectors of the coupling matrix 
M; and the stability of the synchronized states is determined by the corresponding eigenvalues. The set of normal 
modes of the coupling matrix reflect the connectivity of the tree. These modes have complex spatial forms but they 
are analogous to the Fourier eigenmodes arising in regular Euclidean lattices. On the other hand, the distribution of 
eigenvalues of M and their degeneracies are nonuniform. These features affect the bifurcation properties of dynamical 
systems such as coupled map defined on trees. The scaling structure of the synchronized, period-doubled states is 
similar for both uniform and hierarchical lattices, but the nature of the bifurcation boundaries is different. For trees, 
the boundary curves are determined by the spectrum of eigenvalues of the coupling matrix, which has a nonuniform 
density. The nonuniform distribution of eigenvalues leads to gaps or niches in the boundary curves that are not 
present for coupled maps on uniform lattices, where the spectrum of eigenvalues is continuous. 

We have examined only the simplest spatiotemporal patterns that can be formed on treelike geometries; however, the 
formalism presented in Sec. II can be applied to many other processes, such as nontrivial collective behavior, excitation 
waves, phase transitions, domain segregation and growth on trees. The formalism is also useful for continuous-time 
local dynamics. Similarly, extensions of this work arc possible in order to include networks with variable ramification 
and/or depths. 

The study of dynamical systems defined on trees and other nonuniform substrates should allow us to gain in- 
sight into previously unexplored spatiotemporal phenomena on inhomogeneous systems and to understand better the 
relationship between topology and collective properties of networks. 
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fig. 1 



FIG. 1. Tree with ramification R — 3 and depth L — 3, showing the labels on the cells. The corresponding vector-component 
index j is indicated besides each cell. 
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FIG. 2. A tree with R = 3, L = 2. a) The three non-degenerate, layered eigenvectors {u„} and their corresponding 
eigenvalues, b) Degenerate eigenvectors {v^ s } and corresponding eigenvalues. 
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FIG. 3. Spectrum of eigenvalues at increasing depths L, for a tree with ramification R — 3. Eigenvalues {&„} are indicated 
by black dots (•). Other symbols indicate eigenvalues {a ms } as follows: <zi m (A); aim (*); fl3m (x); ct4 m (o); a 5m (+). 
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fig. 4(b) 
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FIG. 4. a) Full spectrum of eigenvalues of the coupling matrix M for a tree with R = 2, L = 11, showing their degeneracy. 
For clarity, eigenvalues {b n } are shown with dots just bellow the zero line. The vertical axis shows the degeneracy of the 
eigenvalues {a ms } divided by N, indicated by vertical bars at each eigenvalue, b) The measure of the set of all eigenvalues of 
M. 



15 




fig. 5(b) 
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FIG. 5. a) The boundary curves S3 = ±1 given by Eqs. (^)-(^) for the period-two, synchronized states of a tree with 
R = 3, L = 3. The curve S3 (64) = +1 is signaled by arrows. The upper curves correspond to the r.h.s equal to —1 for both 
types of eigenvalues. The interior region bounded by these curves is where stable, synchronized, period-two states exist in 
the plane (7, A), b) Magnification of the upper curves in a) showing the gaps in the stability boundary of the period-two, 
synchronized states. Curves corresponding to several eigenvalues are indicated by arrows. The cross just beyond the boundary 
S^(a22) = —1 indicates the parameters 7 and A used in Fig. 6. 




FIG. 6. Inhomogeneous, period-4, state at parameter values 7 = —0.22, A = 3.2367 for a tree with R = 3, L = 3. This 
pattern is a linear combination of the six eigenmodes associated to the eigenvalue 0,22- White corresponds to a zero value and 
black to a value equal to one. 



17 



